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ABSTRACT 


The  transmission  of  a  one-dimensional 
shock  wave  across  an  evacuated  barrier  in  a 
liquid  is  modeled.  In  this  model  the  inter¬ 
action  of  the  shock  wave  with  the  incident 
side  of  the  barrier  causes  the  liquid  to 
cavitate.  The  cavitated  liquid  then  moves 
across  the  barrier  and  collides  with  the 
liquid  on  the  opposite  side  of  the 
barrier ,thus  producing  a  transmitted  wave 
into  the  liquid. 

The  pressure  history  of  the  transmitted 
wave  is  determined  from  that  of  the  incident 
shock  wave.  Although  the  impulse  of  the 
transmitted  wave  is  the  same  as  that  of  the 
incident  shock  wave,  the  transmitted  wave  can 
have  considerably  lower  peak  pressure  and 
energy  flux  density.  Possible  extensions  of 
this  model  are  discussed.  Also,  potential 
military  and  industrial  applications  are 
considered. 


ADMINISTRATIVE  INFORMATION 

This  work  was  funded  by  the  David  Taylor  Research  Center 
Independent  Exploratory  Development  ( IED)  program  under  Work  Unit 
1750-236  for  FY-1990. 


INTRODUCTION 

Surface  ships  and  submarines  must,  be  protected  from  under¬ 
water  explosions  so  that  they  can  survive  to  accomplish  their 
missions.  Furthermore,  it  is  important  that  the  effects  of 
underwater  explosions  in  industrial  applications  (such  as  under¬ 
water  construction  and  demolition)  be  contained  in  the  desired 
region  so  that  machinery  and  personnel  are  protected.  One  method 
of  protection  is  to  place  a  barrier  between  an  explosion  and  an 
object  to  reduce  the  damaging  effects  of  the  shock  wave.  This 
barrier  might  consist  of  a  foam  or  bubble  screen  or  some  more 
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complicated  arrangement.  However,  to  find  the  proper  types  of 
barriers  it  is  first  necessary  to  understand  the  fundamental 
mechanisms  that  reduce  the  damage.  To  achieve  this  it  is  worth¬ 
while  to  examine  and  analyze  simple  models  of  protective  barriers 
to  identify  and  to  model  mathematically  the  physics  of  the  phe¬ 
nomena.  It  is  the  objective  of  this  program  to  develop  such  a 
model.  In  this  model  a  simple  one-dimensional  shock  wave  is 
incident  on  a  barrier  containing  no  material,  i.e.,  a  vacuum. 

The  development  of  the  model  considers  the  cavitation  of  the 
liquid  at  its  interface  with  the  barrier.  This  development  shows 
that  the  decay  in  the  pressure  of  the  incident  wave  causes  a 
reduction  in  the  pressure  and  energy  flux  per  unit  area  of  the 
transmitted  wave.  Of  course,  there  is  no  loss  of  energy  to  the 
vacuum  itself.  Previous  work  on  the  transmission  of  shock  waves 
across  barriers  (composed  of  foam,  air  bubbles,  etc.)  appears  to 
have  concentrated  on  the  transmission  across  the  material  inside 
the  barrier  rather  than  considering  the  effect  of  the  liquid 
cavitating  at  the  interface  of  the  liquid  with  the  barrier,  see, 
e.g.,  Refs.  1-10.  However,  some  work  considering  cavilation  has 
been  performed,  see,  e.g.,  Lyakhov. H  It  is  plausible  that 
existing  experimental  data  that  were  obtained  for  bubble  screens 
can  be  used  to  test  this  theory. 

EVACUATED-BARRIER  CAVITATION  MODEL 
In  the  evacuated-barrier  cavitation  model,  a  one-dimensional  planar 
shock  wave  propagating  in  a  liquid  is  incident  on  a  barrier  that 
is  immersed  in  the  liquid;  see  Fig.  1.  The  pressure  in  the  shock 
wave  is  sufficiently  high  so  that  the  pressure  in  the  barrier  is 
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NOTE  THAT  THESE  ARE  SIMPLY  BOUNDARY  LINES.  THEY 
CONTAIN  NO  MATERIAL  AND  ARE  OF  ZERO  THICKNESS. 


INCIDENT  SHOCK  FRONT 

DIRECTION  OF  PROPAGATION 

- 

EVACUATED 

LIQUID 

BARRIER 

-*g - D - 

Fig.  1.  Evacuated-barrier  model. 
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negligible  in  comparison  to  that  of  the  shock  wave.  Thus,  the 
pressure  within  the  barrier  can  be  considered  as  that  of  a  vacu¬ 
um,  i.e.,  zero.  If  the  incident  pressure  is  sufficiently  high, 
the  reflected  shock  wave  will  causa  a  quasi-continuous  cavitation 
of  the  liquid.  This  will  result  in  a  region  of  cavitation  bub¬ 
bles  between  the  front  of  the  reflected  shock  wave  and  the  barri¬ 
er  as  shown  in  Fig.  2.  The  cavitated  liquid  parcels*  will  then 
move  to  the  opposite  side  of  the  barrier  and  collide  with  it. 
The  shock  wave  resulting  from  rhese  collisions  will  subsequently 
propagate  into  the  liquid  on  the  other  side  of  the  barrier,  as 
shown  in  Fig.  3.  However,  the  times  of  arrival  of  the  liquid 
parcels  are  spread  out  because  of  their  movement  across  the 
barrier.  This  spread  causes  the  shock  wave  that  they  produce  to 
have  a  considerably  lower  pressure  than  that  of  the  incident 
shock  wave.  A  similar  reduction  occurs  in  the  energy  flux  densi¬ 
ty  of  the  shock  wave. 

REVIEW  OF  THEORY  FOR  INITIAL  VELOCITY  OF  CAVITATED  LIQUID 

It  is  instructive  to  review  the  theory  for  determining  the 
velocity  distribution  in  the  liquid  immediate.ly  after  cavitation. 
This  review  is  similar  to  work  that  has  been  reported  by,  e.g., 
Kennard,1^  Cushing, and  Waldo. ^  Consider  a  one-dimensional 
situation  where  a  planar  shock  wave  is  acoustically  propagating 
the  x-  direction  but  in  the  domain  where  x  is  negative  and  is 
incident  on  a  free  surface  at  x  =  0  as  shown  in  Fig.  4.  Also  let 


For  one-dimensional  flow,  these  parcels  are  "slabs,"  as  portrayed  In  Fig.  3  and  for  three-dimensional  flow 
they  might  be  droplets. 
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2. 

Cavitation, 

EVACUATED 

BARRIER 


ORIGINAL  BARRIER  POSITIONS 


MOVING  SLABS  (PARCELS) 


OF  LIQUID 

Fig.  3.  Collision  of  cavitated  slabs  and  transmitted  wave. 
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INCIDENT  SHOCK  FREE  SURFACE 

FRONT  (BARRIER) 


Fig.  4.  Planar  shock  wave  incident  on  surface. 
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the  pressure  of  this  incident  shock  wave  at  the  location  x  and 
time  t  be  given  by 

Pj(x,t)  =  f(ct  -  x),  for  ct  >  x 

=  0,  for  ct  <  x,  (1) 

where  c  is  the  speed  of  sound  in  the  liquid  (see,  e.g.,  Cole,1^ 
p.  143  or  Landau  and  Lifshitz,1®  p.  247).  From  this  equation  it 
is  seen  that  this  shock  wave  arrives  at  the  free  surface  at 
t  =  0.  Now  let  the  pressures  in  the  incident  wave  be  sufficient¬ 
ly  large  so  that  the  pressure  at  the  free  surface  and  the  hydro¬ 
static  pressure  are  negligible  in  comparison.  Immediately  after 
arrival,  a  reflected  wave  will  propagate  from  the  free  surface  in 
opposite  direction  with  the  pressure 

-  f(ct  +  x),  for  ct  £  -x 

and 

0,  for  ct  <  -x; 

see  Fig.  5.  Thus,  if  the  pressures  of  the  incident  wave  are 
positive,  then  the  reflected  wave  is  a  tension  wave.  In  this 
situation,  the  total  pressure  in  the  liquid  will  be  the  pressure 
of  the  incident  wave  plus  the  pressure  of  the  reflected  wave, 
i .e. , 


Pq ( x , t )  =»  f(ct  -  x)  -  f(ct  +  x),  for  ct  >  -  x,  (2) 

where  the  "0"  in  the  subscript  denotes  that  the  reflection  is  at 
the  free  surface  at  x  =  0 .  Note  that  Eq.  2  obeys  the  boundary 
condition,  Pg(0,t)  =  0,  at  the  free  surface  (at  x  =  0).  It  is 
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FREE  SURFACE 
(BARRIER) 


Fig.  5.  Reflected  shock  wave. 
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now  assumed  that  f(z)  is  a  decreasing  function  of  z.  That 
is,  for  any  real  number  z  t  0, 

f(z  +  €)  <  f(z),  for  any  real  number  6  >  0, 
see  Fig.  6.  For  this  decreasing  function  it  is  seen  that  for  x<0 

P0(x,t)  <  0,  for  ct  >  -  x. 

Now  let  the  pressures  in  the  incident  wave  be  sufficiently  large 
so  that  the  magnitude  of  the  (negative)  pressure*  necessary  to 
cause  the  liquid  to  cavitate  is  also  negligible  in  comparison. 
Under  this  condition  the  liquid  will  cavitate  and  form  a  pair  of 
new  free  surfaces  at  x  =  x^  that  are  very  close  to  the  original 
free  surface  (at  x  =  xg  =  0);  see  Eq.  2  for  the  total  pressure  in 
the  liquid.  The  very  small  space  between  this  pair  of  new  free 
surfaces  is  the  first  cavity;  see  Fig.  7.  Defining  p  to  be  the 
mass  density  of  the  liquid  and  using  one-dimensional  linear 
acoustic  theory  (see,  e.g.,  Cole,-^  p.  143  or  Landau  and  Lif- 
shitz,1^  p.  247),  it  is  seen  that  the  velocity  of  the  liquid  at 
the  new  free  surface  that  is  closest  to  the  barrier  on  this 
cavity,  occurring  at  t  =  -x^/c,  is 

V1  *  [ f ( ct  -  x1)  +  f(ct  +  x1)]/(pc) 

-  [f(-2xx)  +  f ( 0 ) ] / (pc ) 

a  2f (-2x1)/(pc) ,  (3) 

because  x^  is  approximately  zero  and  becomes  exactly  zero  in  the 


*'The  problem  for  lower  pressures  is  much  more  complicated,  see,  e.g.,  Prosperetti  [17]. 
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Fig.  6.  Decreasing  incident  shock  wave. 


FIRST 
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PAIR  OF  NEW  FREE 
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BY  CAVITATION 


VELOCITY  AT  NEW 
FREE  SURFACE 


*1 


Fig.  7.  Pair  of  new  free  surfaces  caused  by  cavitation. 
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limit  of  zero  cavitation  pressure  at  the  free  surface.  Because 
the  space  between  the  slabs  is  negligibly  small  immediately  after 
cavitation,  the  location  of  both  surfaces  is  at  x  =  .  It  is 

noted  that,  for  this  first  new  cavitated  slab,  the  approximation 

VX  »  2f ( 0 ) / (pc)  (4) 

could  also  be  used  and  would  be  exact  in  the  limit  of  zero  cavi¬ 
tation  pressure. 

In  similar  fashion,  the  remainder  of  the  incident  wave  will 
reflect  off  of  this  first  cavity  with  the  pressure 

-  f(ct  +  x  -  2xi)  +  Pv,  for  ct  >  -  x  and  x  £  x^ 

where  Pv  is  the  pressure  in  the  cavity  (that  might  be  of  similar 
magnitude  to  the  vapor  pressure  of  the  liquid) .  The  total 
pressure  in  the  liquid  will  be  the  pressure  of  the  incident  wave 
plus  the  pressure  of  the  reflected  wave,  i.e., 

Px(x,t)  =  f ( ct  -  x)  +  .[  —  f ( ct  +  x  -  2xx)  +  Pv], 
for  ct  >  -  x  and  x  £  x^, 

where  the  subscript  "1"  denotes  that  the  reflection  is  at  x  = 

It  is  noted  that  this  equation  obeys  the  boundary  condition, 
Pl(xi,t)  =«  Pv  at  x  =  Now  let  the  pressures  in  the  incident 

wave  be  sufficiently  large  so  that  Pv  is  also  negligible  in 
comparison.  Under  this  condition  the  liquid  will  cavitate  and 
form  a  second  pair  of  new  free  surfaces,  see  Fig.  8,  at  x  =  X2 
that  are  very  close  to  the  first  pair  (at  x  =  x^).  The  veloci¬ 
ty  of  the  liquid  at  the  new  free  surface  that  is  closest  to  the 
barrier  on  this  second  cavity,  occurring  at  t  =  -  X2/C,  is 


SECOND  FIRST 

CAVITY  CAVITY 


x2  X!  Ox 


Fig.  8.  Second  pair  of  new  free  surfaces  caused  by  cavitation. 
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(5) 


V2  =  [f(ct  -  X2)  +  f(ct  +  X2  -  2x1)]/(pc) 

-  [ f ( -2x2 )  +  £(-  2x1)]/(pc) 

»  2f ( -2x2 ) / (pc ) , 

because  X2  =*  and  becomes  exactly  equal  in  the  limit  of  both 
zero  cavitation  pressure  and  Pv. 

As  long  as  the  pressure  in  the  incident  wave  is  sufficiently 
large,  this  process  of  successive  cavitations  and  new  tree  sur¬ 
faces  will  continue  so  that  we  can  consider  the  n-th  cavitation 
in  a  similar  fashion;  see  Fig.  9.  The  incident  wave  will  reflect 
off  of  this  n-th  cavity  with  the  pressure 

-  f(ct  +  x  -  2xn)  +  Pv,  for  ct  >  -x  and  x  <  xn. 

It  is  obvious  from  the  previous  discussion  that  the  total  pres¬ 
sure  obeys  the  boundary  condition  at  x  =  xn.  For  sufficiently 
high  pressures  in  the  incident  wave,  the  liquid  will  cavitate  and 
form  an  n+l-th  pair  of  new  free  surfaces  at  x  =  xn+1  that  is  very 
close  to  the  n-th  pair  (at  x  =  xn) .  The  velocity  of  the  liquid 
at  the  new  free  surface  on  this  n+l-th  cavity  occurring  at 
t  =  -xn+1/c  is 

vn+l  "  tf(ct  -  xn+1)  +  f(ct  +  xn+1  -  2xn)]/(pc) 

“  [f(-2xn+1)  +  f (-  2xn)]/(pc) 

53  2f (-2xn+1)/(pc),  (6) 

because  xn+1  »  xn  and  becomes  exactly  equal  in  the  limit  of  zero 
cavitation  pressure  and  Pv. 

The  liquid  between  these  cavitated  surfaces  can  be  regarded 
as  cavitated  "slabs";  see  Fig.  9.  Because  the  distances  between 
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CAVITATED  SLABS 


Rg.  9.  Pair  of  new  free  surfaces  at  x  =  Xn+1. 
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the  cavitated  surfaces  are  very  small,  the  slabs  are  very 
thin;  and,  because  they  are  very  thin,  the  slabs  have  approxi¬ 
mately  uniform  velocity  throughout  their  thicknesses.  Thus,  the 
approximate  velocity  of  the  n+l-th  slab  given  by  Eq.  6  is  exact 
in  the  limit  where  the  cavitation  pressure  and  Pv  are  zero.  It 
is  mentioned  in  passing  that  Eq.  6  can  also  be  derived  by  divid¬ 
ing  the  impulse  per  unit  area  imparted  by  the  incident  wave  to 
the  n+l-th  slab,  i.e., 

f (“2xn+lH2(xn  -  xn+l)/c]' 
by  the  mass  per  unit  area  of  the  slab,  i.e., 

p(xn  -  xn+1), 

as  was  done  by  Waldo, ^  pp.  2  and  3. 
define  the  velocity  function  V(x)  as 

V(x)  -  2f ( -2x ) / (pc ) , 
from  which 

V(xn)  =*  Vn  (8) 

by  using  Eq.  6  for  the  velocity  of  the  slab. 

TRANSMITTED  WAVE 

After  the  first  slab  cavitates,  it  moves  across  the  barrier 
and  collides  with  the  opposite  side  of  the  barrier;  see  Fig.  10. 
This  collision  produces  an  acoustic  pulse  that  propagates  into 
the  liquid  on  the  opposite  side  of  the  barrier.  In  the  acoustic 
limit,  the  displacement  of  the  liquid  in  the  slab  after  collision 


It  is  now  convenient  to 

(7) 
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is  negligible  so  that  the  position  of  the  end  (i.e.,  the  left- 
hand  side  in  Fig.  10)  of  the  first  slab  remains  at  x  =  D+x-j. 
(remember  that  xj_  <  0)  .  In  like  manner,  after  the  second  slab 
cavitates,  it  also  moves  across  the  barrier  and  collides  with  the 
first  slab  (that  has  previously  collided  with  the  opposite  side 
of  the  barrier  and,  in  the  acoustic  limit,  has  the  same  thickness 
as  when  it  was  formed  by  cavitation);  see  Fig.  11.  This  process 
continues  so  that  after  n-1  succeeding  slabs  pile  up  on  each 
other,  the  n-th  slab  collides  with  the  pile  of  n-1  slabs  at  the 
time, 

(-xn/c)  +  (D/Vn), 

where  D  is  defined  to  be  the  barrier  thickness  (see  discussion  of 
Eq.  6  for  the  velocity  of  the  slab).  The  time  when  the  acoustic 
pulse  from  this  collision  arrives  at  x  =  D  (i.e.,  the  opposite 
side  of  the  original  position  of  the  barrier;  see  Fig.  12)  is 
( see  Eq .  8 ) 

t"n  s  U-Xn/C)  +  (D/Vn)]  +  (-xn/c) 

*  (-2xn/c)  +  [ D/V ( xn ) ] 

*  t'n  +  D/V ( -ct ' n/2 ) ,  (9) 

where 

t'n  “  -2xn/c  (10) 

is  the  arrival  time  of  the  pulse  in  the  limit  of  zero  barrier 
thickness.  It  can  be  seen,  using  Eq.  1  for  the  pressure  of  the 
incident  shock  wave,  that 
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FRONT  OF  COLLISION  PULSE 
FIRST  SLAB  FROM  FIRST  SLAB 

AFTER  COLLISION  / 
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Fig.  10.  Collision  of  first  slab  with  the  opposite  side  of  the  barrier. 


ORIGINAL  SECOND  ‘LENGTH  OF  PULSE 


PULSES  FROM  SLAB  COLLISIONS 
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f(ct'n)  =  p1(0,t'n). 


(11) 


Substituting  Eq.  11  into  Eq.  7,  the  velocity  function,  gives 

V(-ct'n/2)  *  2PI(0,t'n)/(pc).  (12) 

Substituting  Eq.  12  into  Eq.  9,  for  the  time  when  the  acoustic 
pulse  from  this  collision  arrives  at  x  =  D,  yields 

t"n  3  t'n  +  %pcD/PI( 0,t'n).  (13) 

From  this  equation,  it  is  seen  that  it  is  appropriate  to  define 
the  functional  t''  such  that,  in  the  limit  of  zero  slab  thick¬ 
ness  , 


t"n  =  t"(t'n),  (14) 

or  simply 

t"  =  t"(t').  (15) 

Differentiation  of  Eq.  15  yields 

dt'  1  =  1  -  ijpcD[PI(0,t'  )]-2  d _ 1PI(0,t').  (16) 

dt '  dt ' 


This  equation  implies  that  dt''/dt'  >  1,  for  D  >  0,  because 
Pj(0,t' )  is  a  decreasing  function  of  t' .  This  result  is  expected 
because  the  fact  that  the  pressure  of  the  incident  shock  wave  de¬ 
creases  with  time  implies  that  the  cavitated  slabs  have  succes¬ 
sively  decreasing  velocities.  This  causes  a  spread  in  their 
arrival  times  at  the  opposite  side  of  the  barrier  resulting  in 
dt' ' /dt'  >  1 . 


20 


If  the  slabs  were  of  finite  thickness,  then  the  transmitted 
wave  would  consist  of  a  series  of  pulses,  each  corresponding  to  a 
slab  collision;  see  Fig.  12.  However,  in  the  limit  of  zero 
thickness  slabs  (that  occur  for  sufficiently  high  pressure  inci¬ 
dent  shock  waves,  as  explained  previously)  these  pulses  are  of 
zero  duration  and  of  infinite  frequency.  In  this  situation  the 
pulses  are  rapidly  damped  out  by  viscous  effects  so  that  the 
transmitted  wave  becomes  a  smooth  function  of  time  such  that  the 
momentum  of  each  pulse  is  conserved,  but  there  might  be  a  signif¬ 
icant  loss  of  energy.  It  is  emphasized  that  the  loss  of  energy 
occurs  in  the  liquid,  as  it  must  because  there  is  a  vacuum  within 
the  barrier  itself.  Using  Eq.  10  for  the  pulse  arrival  time  in 
the  limit  of  zero  barrier  thickness  and  Eq.  12  for  the  velocity 
function,  it  can  be  easily  seen  that  the  momentum  per  unit  area 
of  the  n+l-th  slab  is 


p(xn  -  xn+1)Vn+1  *  Pi(0,t'n+1)(t'n+1  -  t'n) 

=  Pj(0,t')dt',  (17) 

where  t'  *  t'n+1  an<*  =  t'n+1  “  t'n  -*-n  t^ie  of  zero 

thickness  slabs.  Integration  of  this  equation  over  time  to  a 
given  value  of  t'  provides  the  total  momentum  of  the  collided 
slabs  up  to  t' .  This  total  momentum  is  equal  to  the  total  im¬ 
pulse  of  the  incident  shock  wave  up  to  t'.  This  total  momentum 
must  be  equal  to  the  total  impulse  up  to  t''  =»  t''(t')  of  the 
transmitted  wave  at  x  =  D,  i.e.. 


0,t)dt 


pt'  ' 

PT(D,t)dt, 

0 
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where  PT(x,t)  is  the  pressure  of  the  transmitted  wave  at  any 
location,  x,  and  at  any  time,  t;  see  Eq.  15.  For  this  equality 
of  momentum  and  impulse  to  be  valid  for  any  value  of  t',  it  is 
required  that 

PI(0,t')dt'  »  PT(D,t' ' )dt" 
or 

PT(D,t")  =  PI(0,t'  )/(dt"/dt' ) ,  (18) 

where  dt''/dt'  is  given  by  Eq.  16.  It  is  noted  that  PT(D,t'')  is 
assumed  to  vary  smoothly  from  t''n+1  to  t' 'n+2  rather  than  with 
the  short  pulse  from  t''n+1  to  t '  '  n+1  +  2 ( xn-xn+1 ) /c  that  would 
occur  if  there  were  no  high  frequency  damping  due  to  viscosity; 
see  Fig.  13  and,  e.g.,  Landau  and  Lifshitz,1^  p.  300.  For  a 
given  value  of  t',  the  corresponding  value  of  t' '  now  can  be 
calculated  using  Eq.  13  and  the  transmitted  pressure  at  time  t' ' 
then  can  be  calculated  using  Eq.  18.  Equation  18  shows  that  the 
pressure  of  the  transmitted  wave  at  time  t''  is  reduced  from  that 
of  the  incident  wave  at  the  corresponding  time  t'  by 

(dt' '/dt' J*1. 

For  D  =  0,  there  would  be  no  reduction  in  the  pressure  of  the 
incident  wave  because  dt''/dt  =  1.  However,  for  D  >  0,  there 
would  be  a  reduction  in  the  pressure  of  the  incident  wave  because 
dt''/dt'  >  1,  as  previously  discussed.  It  is  emphasized  that 
there  would  be  no  reduction  in  the  incident  wave  if  the  incident 
wave  were  constant  in  time  because  dt''/dt'  *  1  in  this  case. 

i 

This  accentuates  the  fact  that  in  this  model  all  the  reductions 
in  the  pressures  of  the  incident  shock  wave  are  the  result  of  it 
being  a  decreasing  function  of  time. 
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For  a  three-dimensional  situation,  the  cavitated  slabs  are 
very  unstable  and  would  rapidly  break  up  into  droplets  that 
would  collide  with  the  opposite  side  of  the  barrier.  This  would 
not  result  in  much  change  in  the  velocities  of  the  liquid  in  the 
slabs,  but  it  would  produce  many  more  collisions  and  acoustic 
pulses  that  would  be  even  more  rapidly  damped  than  in  the  one¬ 
dimensional  case.  This  phenomenon  would  produce  a  smoother 
transmitted  wave  for  a  given  slab  thickness  than  that  for  the 
one-dimensional  situation  presented  above. 

ENERGY  TRANSMISSION  RATIO 

The  energy  flux  per  unit  area  (i.e.,  the  energy  flux  densi¬ 
ty)  of  the  incident  wave  is  determined  using 

Ex  =  (pc)-lffp1(0,t')]2dt',  (19) 

J0 

see  Cole,-1-^  p.  14  3.  In  like  manner,  the  energy  flux  per  unit 
area  of  the  transmitted  wave  is  determined  using 

Et  =  (pc)“1jfpT(D,t")]2dt".  (20) 

Substituting  Eq.  18  for  PT  in  Eq.  20,  and  using 
dt"  =  (dt ' '  /dt'  )dt ' ,  gives 

Et  -  (pc)-1 ffpI(0,t' )  ]2(dt"  /dt' )-1dt' .  (21) 

J0 

Equation  21  shows  that  the  energy  flux  density  of  the  transmitted 
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wave  is  less  than  that  of  the  incident  wave  because  dt''/dt'  >  1; 
see  the  discussion  of  Eq.  16.  The  energy  transmission  ratio, 
Et/Ej,  now  can  be  determined  by  using  Eqs .  19  and  21.  As  ex¬ 
plained  above,  the  impulse  of  the  transmitted  wave  is  the  same  as 
that  of  the  incident  wave. 

SPECIAL  CASES 

To  demonstrate  the  use  of  this  model,  the  special  cases  of 
the  exponentially  and  hyperbolically  decaying  incident  shock 
waves  are  now  presented. 

Case  1 .  Exponentially  Decaying  Incident  Shock  Wave 

The  pressure  for  the  special  case  of  an  exponentially  decay¬ 
ing  incident  shock  wave  is 

PI(0,t')  =  PIQexp(-t' /9 ) ,  (22) 

where  Pjg  is  the  peak  pressure  (occurring  at  t'  =0)  and  9  is  the 
exponential  time  decay  constant.  In  this  case,  Eq.  16  becomes 

dt> '  =  1  +  %pcD  exp(t'/9)/(PI09) 
dt' 

-  1  +  %D*exp(t'/0) ,  (23) 


where 


D*  =  pcD/(PlO0)  (24) 

is  the  nondimensional  barrier  thickness .  For  a  given  value  of 
t' ,  the  corresponding  value  of  t''  now  can  be  calculated  using 
Eq .  13,  and  the  transmitted  pressure  at  time  t' ’  then  can  be 
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calculated  using  Eqs.  18  and  23.  For  t'  =0 

t"  =  %pcD/Pl0  =  ^D*9  (25) 

and  the  initial  (peak)  pressure  of  the  transmitted  wave  is 

PT(D,t")  =  PT(D,%D*9)  =  PI0/(1  +  %D*).  (26) 

Equation  26  shows  that  the  peak  pressure  of  the  transmitted  wave 
is  less  than  that  of  the  incident  wave  by  a  factor  of  (l+^D*)-1. 
For  sufficiently  large  values  of  t',  it  can  easily  be  seen  from 
Eqs.  13  (for  t''),  22  (for  Pj),  and  24  (for  D*)  that 

t"  *  *iD*9  exp ( t '  / 8  )  (27) 

and,  from  Eqs.  18  and  23,  that 

PT(D,t")  **  2PJQD*"1  exp( -2t ' / 0 ) .  (28) 

For  sufficiently  large  values  of  t'',  Eqs.  27  and  28  imply  that 

PT(D,t"  )  »  i*D*(9/t"  )2PI0.  (29) 

Equation  29  shows  that  for  large  times  the  transmitted  wave 
decays  as  the  reciprocal  of  the  square  of  the  time.  This  is  a 
much  slower  decay  than  that  of  the  incident  wave  that  decays 
exponentially  with  time. 

Fro;  i  Eq.  19,  the  energy  flux  per  unit  area  of  the  incident 
shock  wave  becomes 

Ej  =  (pc)'1  Pi029/2,  (30) 
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and,  from  Eq.  21,  the  energy  flux  per  unit  area  of  the  transmit¬ 
ted  wave  becomes 

Et  =(pc)_1 r“plO2exp(-2t'/0) [l+%D*exp(t'/9) ]-1dt' ,  (31) 

J  0 

by  using  Eq.  23.  Substituting  t'  =  9  logeu  into  Eq.  31  and 
dividing  the  result  by  Eq.  30,  the  energy  transmission  ratio  for 
this  case  can  be  expressed  as 

Et/Ez  =  r°°2  (u3  +  ^3D*u4)"1du.  (32) 

J1 

From  this  equation,  in  the  extreme  situation  where  D*  *  0, 

Et/Et  =  1,  (33) 

as  expected.  Also,  in  the  extreme  situation  where  D*  becomes 
arbitrarily  large,  it  can  be  seen  from  Eq.  32  that 

Ex/Et  =  4/ (3D*).  (34) 

It  can  be  easily  verified  for  these  two  extreme  situations  (i.e., 
Eqs.  33  and  34)  that  the  expression 

(1  +  3D*/4)-1 

has  the  same  asymptotic  behavior.  Also,  by  numerically  integrat¬ 
ing  Eq.  32,  it  can  be  shown  that  for  0  <  D* 

(1  +  3D*/4)_1  -  5  <  Et/Ej  <  (1  +  3D*/ 4 ) -1  , 
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where*  6  *  0.0261.  Thus,  because  of  the  correct  asymptotic 
behavior  and  the  small  value  of  6,  it  appears  that  for  many 
purposes  the  approximation, 

Et/Ex  -  (1  +  3D*/4)-1,  (35) 

should  be  adequate  for  this  case. 

Thus,  for  the  exponentially  decaying  incident  shock  wave, 
the  peak  pressure  and  energy  flux  density  of  the  transmitted  wave 
are  reduced,  but  the  transmitted  wave  does  not  decay  as  rapidly 
as  does  the  incident  wave.  This  is  expected  because  the  impulse 
of  the  transmitted  wave  is  always  equal  to  that  of  the  incident 
wave  for  this  model. 

Case  2.  Hvoerbolicallv  Decaying  Incident  Shock  Wave 

In  the  special  case  of  a  hyperbolically  decaying  incident 
shock  wave, 

Pj^t')  -  PI0(1  +  t'/Q)”1,  for  0  S  t'  <  x  (36) 

=  0,  otherwise 

where,  as  in  the  previous  special  case,  Pjq  is  the  peak  pressure 
(occurring  at  t'  =0),  0  is  the  hyperbolic  time  decay  constant 
(defined  in  analogous  fashion  to  the  previous  special  case)  and  x 
is  the  cut-off  time.  In  this  case,  Eq.  16  becomes 

dt>  *  =  1  +  JjD*,  for  0  £  t'  <  x  (37) 

dt' 


*'Note  that  at  D*  -  0.48,  the  numerical  Integration  of  Eq.  32  yields  ET/EI  •  0.7092.  This  Is  approximately 
equal  to  the  lower  bound  (l.e.,  the  left-hand  side)  of  this  Inequality. 
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where  D*  is  defined  to  be  the  same  as  in  Eq.  24.  However,  unlike 
the  previous  special  case  (Eq.  23),  Eq.  37  is  constant  in  time 
except  for  the  discontinuity  at  t'  =  x.  For  this  hyperbolically 
decaying  shock  wave,  Eq.  13  becomes 


t"  =  ijD*0  +  (1  +  *sD*)t',  for  0  S  t'  <  x 

for  t'  >  r 


=  oo  , 


(38) 


and  Eq.  18  becomes 

PT(D,t")  =  PI0(1  +  t,/0)“1/(l  +  JsD*),  for  0  <  t'  <  x 

=  0,  otherwise.  (38a) 


Solving  for  t'  in  Eq.  38  and  substituting  the  resulting  expres¬ 
sion  in  Eq.  38a  gives 


PT(D,t' ' ) 


PI0  t"-JjD*0 

-  [  1  +  - 

1  +  *5D*  8(l+JjD*) 


]" 


1 


=  PI0(1  +  ^D*  +  (t"  -  hD *0)/0]"1 

=  PI0(1  +  t"/0)_1,  for  0  £  t ' '  — %D*0  <  x ( 1+%D* ) 
=  0,  otherwise.  (39) 

For  this  case  Eq.  39  shows  that  when  0  5  t''-^D*0  <  x(l+35D*),  the 
transmitted  wave  has  the  same  form  as  the  incident  wave ;  see  Eq. 
36.  However,  because  the  pressure  of  the  transmitted  wave  is 
zero  for  t''  <  %D*0,  the  peak  pressure  of  the  transmitted  wave  is 
decreased.  Also,  the  duration  in  increased.  The  transmitted  wave 
is  attenuated  in  this  fashion. 
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For  the  hyperbolical ly  decreasing  shock  wave,  Eq.  19  becomes 


Ex  =  (pc)"1  PI029[1  -  (1  +  x/9 ) '”1  ] ,  (40) 

Eq.  21  becomes 

ET  =  (pc)"1  PIO20[1  -  (1  +  x/8)“1]/(l  +  JjD*),  (41) 

and  the  energy  transmission  ratio  becomes 

ET/EX  =  (1  +  JsD*)-1.  (42) 

Note  that  this  ratio  is  independent  of  x.  Comparing  Eq.  42  with 
Eq.  35  shows  that  the  energy  transmission  ratio  for  the  hyperbol- 
ically  decreasing  shock  wave  is  greater  than  that  of  the  exponen¬ 
tially  decreasing  shock  wave. 

LOWER  PRESSURE  INCIDENT  SHOCK  WAVES 

The  theoretical  development  that  was  presented  above  can  be 
easily  extended  to  the  problem  nf  ir^i^er.t  check  waves  with 
lower  pressures  and  barriers  with  finite  pressures.  This  can  be 
accomplished  with  simple  modifications  of  Eqs .  1-8  to  treat 
cavitated  slabs  with  finite  thicknesses,  assuming  that  the  liquid 
cavitates  at  a  definite  pressure.  However,  the  situation  could 
be  much  more  complicated  because  one  must  explicitly  consider  the 
dynamics  of  the  growth  of  the  cavitation  bubbles  and  the  effect 
of  impurities  (and  their  statistical  distributions)  and  other 
phenomena;  see,  e.g.,  Prosperetti . 17  Also,  the  treatment  of  the 
transmitted  pulse  can  be  easily  modified  in  a  similar  fashion. 
Of  course,  for  sufficiently  low  pressure  shock  waves  that  are 
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incident  on  a  barrier  that  contains  material,  the  usual  acoustic 
treatment  of  sound  wave  reflection  and  transmission  can  be  ap¬ 
plied;  see,  e.g.,  Landau  and  Lifshitz,1^  p.  253.  However,  these 
modifications  require  more  input  parameters  and  might  not  be  of 
much  practical  use  because  protection  is  most  important  for 
incident  shock  waves  with  very  high  pressures.  Although  this 
phenomenon  should  scale  for  sufficiently  high  pressures,  it  might 
not  scale  for  low  pressures.  This  is  because  the  effect  of 
surface  tension  and  viscosity  might  be  very  important  during 
cavitation  and  collision  for  sufficiently  low  pressures. 

EXTENSION  OF  MODEL  TO  THREE  DIMENSIONS 

Perhaps  the  most  elementary  three-dimensional  extension  of 
the  evacuated-barrier  cavitation  model  is  a  model  of  an  incident 
shock  wave  emanating  from  an  explosion  that  can  be  considered  as 
a  point  source.  It  is  noted  that  this  situation  might  only  be 
two-dimensional  because  of  possible  axial  symmetry  about  the  line 
from  the  point  source  to  the  nearest  point  on  the  barrier.  For 
this  situation,  a  simple  model  can  be  applied  for  the  initial 
velocity  distribution  after  the  cavitation  (as  shown  in  Fig.  9) 
for  the  interaction  of  an  incident  shock  wave  with  an  arbitrary 
angle  of  incidence;  see,  e.g.,  Waldo14  and  cf.  Eq.  7. 

Using  this  velocity  distribution,  the  pressure  history  on 
the  opposite  side  of  the  barrier  can  be  determined  by  employing 
the  same  concepts  that  were  developed  previously  for  the  one¬ 
dimensional  situation  in  the  evacuated-barrier  cavitation  model. 
The  calculation  of  the  pressure  history  for  any  distance  beyond 
the  opposite  side  of  the  barrier  would  require  a  solution  of  the 
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acoustic  wave  equation.  The  boundary  condition  for  this  solution 
is  the  pressure  history  at  each  location  on  the  opposite  side  of 
the  barrier  (as  calculated  using  the  model  for  arbitrary  angles 
of  incidence  that  was  described  previously) .  However,  even 
without  explicitly  solving  this  problem,  it  appears  to  be  reason¬ 
able  that  the  pressures  on  the  opposite  side  of  the  barrier  would 
rapidly  decrease  for  locations  with  increasing  angle  of  inci¬ 
dence.  This  would  mean  that  the  pressures  of  the  transmitted 
wave  would  rapidly  decrease  with  increasing  distances  beyond  the 
opposite  side  of  the  barrier  (even  along  the  line  of  normal 
incidence  of  the  shock  wave,  i.e.,  the  axis  of  symmetry).  Thus, 
it  can  be  speculated  that  there  would  be  a  decrease  in  the  im¬ 
pulse  of  the  transmitted  shock  wave  at  a  given  location  from  what 
it  would  have  been  at  that  location  if  there  were  no  barrier. 
However,  a  full  treatment  of  the  problem  would  be  necessary  to 
justify  this  speculation. 

For  explosions  that  are  not  close  to  the  barrier,  the 
bubble  that  is  produced  by  the  explosion  will  have  most  of  its 
expansion  take  place  considerably  after  the  shock  wave  interac¬ 
tion.  Cole,1^  p.  312,  has  shown  that  a  bubble  will  be  repelled 
by  a  free  surface  such  as  that  at  the  evacuated  barrier  and  will 
be  attracted  by  an  rigid  surface.  Thus,  it  is  plausible  that  the 
bubble  repulsion  by  the  evacuated  barrier  will  reduce  the  effect 
of  the  shock  wave  that  emanates  from  the  collapse  of  the  bubble 
simply  because  the  collapse  occurs  further  away  than  if  the 
evacuated  barrier  were  not  present.  Also  the  pressure  of  the 
transmitted  3hock  wave  produced  by  the  collapse  of  the  bubble  is 
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reduced  by  the  same  effect  as  was  modeled  in  the  development 
previously,  see  Eqs .  13,  16  and  18.  Unfortunately,  a  circum¬ 
stance  might  occur  in  which  the  explosion  is  too  close  to  the 
barrier,  the  barrier  is  not  sufficiently  thick,  and  the  object  to 
be  protected  by  the  barrier  is  too  close  to  the  opposite  side  of 
the  barrier.  In  this  circumstance,  the  bubble  expansion  might  be 
sufficiently  large  to  penetrate  the  barrier  and  to  interact  with 
the  object.  If  the  object  is  solid  and  cannot  have  much  move¬ 
ment,  the  bubble  could  be  attracted  to  the  barrier  and  collapse 
on  the  object.  This  collapse  could  cause  unacceptable  damage  to 
the  object.  However,  even  in  this  circumstance  the  damage  could 
be  substantially  reduced  by  the  presence  of  the  barrier. 

EXTENSION  OF  MODEL  TO  BARRIERS  CONTAINING  MATERIAL 
Consider  a  barrier  containing  a  material  of  very  low  density 
and  very  low  crushing-resistance  pressure.  It  is  plausible  that 
the  reductions  in  the  pressures  of  the  transmitted  shock  waves 
caused  by  this  barrier  could  be  about  the  same  as  those  for  an 
evacuated  barrier.  However,  the  thickness  of  this  evacuated 
barrier  would  be  less  than  the  actual  thickness  of  the  barrier 
containing  the  material.  This  "effective  thickness"  might  be 
equal  to  that  of  the  vacant  space  remaining  in  the  barrier  if 
the  material  were  compressed  so  that  it  would  have  a  mass  density 
that  is  equal  to  that  of  the  liquid.  Such  materials  might  in¬ 
clude  gasses  (e.g.,  air,  carbon  dioxide,  water  vapor  or  steam, 
etc.),  foams  (see,  e.g.,  Kudinov  et  al .  ^ and  Lyakhov,^-1), 


bubble  screens,  styrofoam,  solid  foams  made  of  elaborate  compos¬ 
ite  materials,  some  complicated  combination  of  materials,  or 
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layers  of  materials.  However/  it  would  not  be  surprising  if  the 
reductions  of  pressures  in  the  transmitted  shock  waves  were 
greater  than  those  predicted  by  an  evacuated  barrier  with  an 
effective  thickness  obtained  in  this  manner.  The  greater  reduc¬ 
tions  might  occur  because  of  the  (possible)  absorption  of  the 
incident  shock  wave  energy  by  the  material  in  the  barrier.  Also 
it  is  possible  that  the  barrier  could  be  behaving  as  if  it  were 
evacuated  throughout  its  full  unadjusted  thickness,  without  the 
adjustment  of  the  thickness  as  described.  Of  course,  it  is  possi¬ 
ble  that  a  barrier  could  allow  more  energy  transmission  than  that 
of  the  evacuated-barrier  cavitation  model.  More  energy  transmis¬ 
sion  would  occur  in  the  trivial  case  where  the  barrier  is  filled 
with  the  same  liquid  as  that  outside  the  barrier.  In  this  case 
there  would  be  no  change  in  the  transmitted  shock  wave  from  that 
of  the  incident  shock  wave. 

More  dense  materials  (e.g.,  steel,  concrete,  sand,  etc.) 
might  reflect  a  substantial  amount  of  the  energy  of  the  incident 
shock  wave  and  reduce  the  pressures  of  the  transmitted  shock 
wave.  But  barriers  composed  of  such  dense  materials  are  fre¬ 
quently  impractical.  Of  course,  more  elaborate  mathematical 
models  of  the  interaction  of  shock  waves  with  barriers  composed 
of  various  materials  can  be  formulated  using  the  concepts  consid¬ 
ered  in  the  evacuated-barrier  cavitation  model.  Also,  further 
considerations  such  as  the  drag  forces  on  the  liquid  droplets  as 
they  move  across  the  barrier  might  be  included.  A  simple  barrier 
with  negligible  mass  density  that  uniformly  crushes  with  a  con¬ 
stant  pressure  can  be  easily  treated  using  the  concepts  in  Ref . 14 . 
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EXPERIMENTS  TO  TEST  THE  THEORY 


To  test  this  theory,  see  Eqs .  13,  16,  and  18,  experiments 
could  be  performed  for  a  one-dimensional  evacuated  barrier  by 
using  an  underwater  shock  tube.  Also,  calculations  using  the 
theory  could  be  compared  with  existing  data  for  bubble  screens 
and  low-density  foams.  It  would  not  be  expected  that  these 
calculations  would  agree  precisely  with  the  data.  But  it  is 
plausible  that  the  evacuated-barrier  model,  with  an  effective 
thickness  (as  previously  discussed),  might  provide  a  good  first 
approximation . 


POSSIBLE  APPLICATIONS 

Because  almost  all  damaging  shock  waves  decay  with  time 
(Cole,15  p.  110),  the  evacuated-barrier  concept  could  be  used  to 
protect  surface  ships  and  submarines  by  encasing  the  hull  with  a 
very  low  density  solid  foam  (as  previously  discussed).  However, 
this  foam  would  have  to  be  strong  enough  to  withstand  conditions 
at  sea.  In  addition,  for  submarines  the  foam  would  have  to 
resist  the  pressures  at  the  collapse  depth  of  their  pressure 
hulls.  Such  foams  might  add  too  much  to  the  size  and  weight  of 
the  vessel  to  be  practical. 

It  is  conceivable  that  a  bubble  screen  could  be  employed  for 
naval  ships  and  submarines  during  an  attack.  However,  the  prob¬ 
lems  of  additional  piping  and  supply  of  the  gas  for  the  bubbles 
would  have  to  be  surmounted.  Steam  from  the  boilers  could  be 
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used  if  enough  could  be  generated  without  severely  hindering  the 
mobility  of  the  craft.  Also,  a  chemical  could  be  released  by  the 
craft  that  would  react  with  the  water  and  form  bubbles  that  could 
provide  an  effective  screen.  This  chemical  could  be  injected 
into  the  water  as  coated  pellets  at  the  bow  of  the  craft  when 
under  attack.  The  coatings  would  dissolve  at  different  times  and 
thus  the  chemical  would  react  at  different  times.  This  would 
provide  a  bubble  screen  for  the  full  length  of  the  hull. 

There  are  many  other  schemes  that  could  be  investigated,  and 
numerous  industrial  applications  of  bubble  screens  such  as  the 
protection  of  underwater  structures  from  explosions  during  con¬ 
struction  and  demolition. 

Note  that  because  the  impulse  is  not  reduced  by  the  barrier 
(see  the  discussion  preceding  Eq.  18),  there  might  not  be  much 
reduction  in  damage  imparted  to  equipment  that  responds  mainly  to 
impulse  rather  than  peak  pressure  or  energy  flux  density.  As 
previously  discussed,  the  bubble  produced  by  the  explosion  might 
collapse  on  the  barrier,  penetrate  it,  and  cause  additional  dam¬ 
age.  However,  because  of  the  low  density  of  the  barrier,  the 
bubble  might  be  repelled  by  it  and  not  penetrate  the  barrier; 
see,  e.g.,  Cole,1^  pp.  312-352. 

SUMMARY 

The  evacuated-barrier  cavitation  model  was  developed.  In 
this  model  a  one-dimensional  shock  wave  in  a  liquid  is  incident 
on  the  surface  of  an  evacuated  barrier.  Because  of  this  inter¬ 
action,  the  liquid  cavitates.  The  slabs  of  cavitated  liquid 
subsequently  move  across  the  barrier  and  collide  with  the  liquid 
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surface  on  the  opposite  side  of  the  barrier.  These  collisions 
produce  a  high  frequency  of  pulses  in  the  liquid  on  the  opposite 
side  of  the  barrier  that  are  rapidly  damped  out.  The  pressures 
and  energy  flux  density  of  the  resulting  wave  that  is  transmitted 
into  the  liquid  on  the  opposite  side  of  the  barrier  can  be  con¬ 
siderably  reduced.  These  reductions  are  the  result  of  the  fact 
that  the  decay  in  the  pressure  of  the  incident  shock  wave  causes 
the  distances  between  the  cavitated  slabs  to  spread  as  they  move 
across  the  barrier.  Spreading  causes  the  transmitted  wave  to 
have  lower  pressures  and  longer  duration  than  the  incident  shock 
wave.  The  impulse  of  the  transmitted  wave  is  the  same  as  that  of 
the  incident  wave..  The  pressure  history  of  the  transmitted  wave 
can  be  determined  for  a  decreasing  incident  wave  using  Eqs.  13/ 
16,  and  18.  For  the  important  case  of  an  exponentially  decaying 
shock  wave,  the  peak  pressure  of  the  transmitted  wave  is  given  by 
Eq.  26  and  the  energy  flux  density  of  the  transmitted  wave  can  be 
determined  using  Eq.  32  or  its  approximate  form,  Eq.  35.  For  a 
hyperbolically  decaying  incident  shock  wave  it  was  shown  (Eq.  39) 
that  the  transmitted  wave  function  has  the  same  form  as  that  of 
the  incident  wave  but  is  attenuated.  Extensions  of  the 
evacuated-barrier  cavitation  model  to  treat  three-dimensional 
situations  and  to  include  various  materials  within  the  barrier 
were  discussed.  Applications  of  the  model  for  the  shock  protec¬ 
tion  of  surface  ships  and  submarines  as  well  as  industrial  appli¬ 
cations  were  considered. 
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